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Abstract 

Given a hypergraph H with m hyperedges and a set A of m pins, 
i.e. globally fixed subspaces in Euclidean space R'*, a pinned subspace- 
incidence system is the pair {H,X), with the constraint that each pin in 
X lies on the subspace spanned by the point realizations in R"^ of vertices 
of the corresponding hyperedge of H. We are interested in combinatorial 
characterization of pinned subspace-incidence systems that are minimally 
rigid, i.e. those systems that are guaranteed to generically yield a locally 
unique realization. As is customary, this is accompanied by a characteri¬ 
zation of generic independence as well as rigidity. 

In a previous paper m, we used pinned subspace-incidence systems 
towards solving the fitted dictionary learning problem, i.e. dictionary learn¬ 
ing with specified underlying hypergraph, and gave a combinatorial char¬ 
acterization of minimal rigidity for a more restricted version of pinned 
subspace-incidence system, with H being a uniform hypergraph and pins 
in X being 1-dimension subspaces. Moreover in a recent paper [2], the spe¬ 
cial case of pinned line incidence systems was used to model biomaterials 
such as cellulose and collagen fibrils in cell walls. In this paper, we extend 
the combinatorial characterization to general pinned subspace-incidence 
systems, with H being a non-uniform hypergraph and pins in X being 
subspaces with arbitrary dimension. As there are generally many data 
points per subspace in a dictionary learning problem, which can only be 
modeled with pins of dimension larger than 1, such an extension enables 
application to a much larger class of fitted dictionary learning problems. 


1 Introduction 

A pinned subspace-incidence system {H,X) is an incidence constraint system 
specified as a hypergraph H together with a set X of subspaces or pins in 
in one-to-one correspondence with the hyperedges of H. A realization of 
{H,X) assigns points in to the vertices of H, thereby subspaces to the hy¬ 
peredges of H. The subspace in corresponding to a hyperedge of H contains 
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the associated pin from X. We are interested in characterization of pinned 
subspace-incidence systems that are minimally rigid, i.e. those systems that are 
guaranteed to generically yield a locally unique realization. 

In a previous paper m, we used pinned subspace-incidence systems towards 
solving/itted dictionary learning problems, i.e. dictionary learning with specified 
underlying hypergraphs. Dictionary learning (aka sparse coding) is the problem 
of obtaining a set dictionary vectors that sparsely represent a set of given data 
points in Geometrically, such a sparse representation can be viewed as a 
subspace arrangement spanned by the dictionary vectors that contains all the 
data points. In fitted dictionary learning, the underlying hypergraph H of the 
subspace arrangement is specified, and the problem becomes a pinned subspace- 
incidence systems with the pins corresponding to the span of data points on each 
subspace. 

Moreover in a recent paper [5], we have used pinned subspace-incidence 
systems in modeling biomaterials such as cross-linking cellulose and collagen 
microfibrils in cell walls 0II m. In such materials, each fibril is attached 
to some fixed larger organelle/membrane at one site, and cross-linked at two 
locations with other fibrils. Consequently, they can be modeled using a pinned 
line-incidence system with H being a graph, where each fibril is modeled as an 
edge of H with the two cross-linkings as its two vertices, and the attachment is 
modeled as the corresponding pin. 

We gave in m a combinatorial characterization of minimal rigidity for a 
restricted version of pinned subspace-incidence system, with the underlying 
hypergraph H being a uniform hypergraph and pins in X being 1-dimension 
subspaces. 

2 Contributions 

In this paper, we extend the combinatorial characterization of minimal rigidity 
to general pinned subspace-incidence systems, where H can be any non-uniform 
hypergraph, and each pin in X is a subspace with arbitrary dimension. Such an 
extension enables application to a much larger class of fitted dictionary learning 
problems, since there are generally many data points per subspace in a dictio¬ 
nary learning problem, which can only be modeled with pins of dimension larger 
than 1. 

As in our previous paper m, we apply the classic method of White and 
Whiteley m to combinatorially characterize the rigidity of general pinned 
subspace-incidence systems. The primary technique is using the Laplace decom¬ 
position of the rigidity matrix, which corresponds to a map-decomposition [m 
of the underlying hypergraph. The polynomial resulting from the Laplace de¬ 
composition is called the pure condition, which characterizes the conditions that 
the framework has to avoid for the combinatorial characterization to hold. 

Previous works on related types of geometric constraint frameworks include 
pin-collinear body-pin frameworks [9], direction networks |20) . slider-pinning 
rigidity [16] , body-cad constraint system [8] , Ic-frames [HI [19] , and affine rigid- 
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(a) 



Figure 1: (a) A rank-3 non-uniform hypergraph with 6 vertices and 5 edges, 
where si = S 2 = 3, S 3 = S 4 = 2, S 5 = 1. (b) A pinned subspace-incidence 
framework in d = 4 (projectivized in P^(R)) on the hypergraph from (a), with 
mi = m 3 = 2, TO 2 = 3, 7714 = 7715 = 1, where the crosses on each hyperedge 
represent points spanning the associated pin. 


ity m- However, we are not aware of any previous results on systems that are 
similar to pinned subspace-incidence systems. 


3 Preliminaries 

In this section, we introduce the formal definition of pinned subspace-incidence 
systems and basic concepts in combinatorial rigidity. 

A hypergraph H = (H, E) is a set V of vertices and a set E of hyperedges, 
where each hyperedge is a subset of V. The rank r{H) of a hypergraph H is 
the maximum cardinality of any edge in E, i.e. r[E[) = maxg^^gs s(efe), where 
s(efc) denotes the cardinality of the hyperedge e^. A hypergraph is s-uniform 
if all edges in E have the same cardinality s. A configuration or realization of 
a hypergraph H = {V,E) in is a mapping of points in to the vertices of 
E[, p : V ^ When there is no ambiguity, we simply use pi to denote the 
point p(vi), p{ek) to denote the set of points {p{vi)\vi G Ck}, and Sk to denote 
the cardinality s(efc). 

An example of a rank-3 hypergraph is given in Figure fTal 

In the following, we use (S) to denote the subspace spanned by a set S of 
points in R'^. 

Definition 1 (Pinned Subspace-Incidence System). A pinned sub space-ineidence 
system in R"^ is a pair {H,X), where E[ = {V,E,m) is a weighted hypergraph 
with hyperedges of rank r(H) < d, and X = {a;i,a; 2 ,... ,x\e\} is a set of pins 
(subspaces of R*^) in one-to-one correspondence with the hyperedges of H. Here 
the weight assignment is a function m : E{H) Z+, where m(e) denotes the 
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dimension of the pin associated with the hyperedge e. Often we ignore the 
weight m and just refer to the hypergraph {V,E) as H. 

A pinned subspace-incidence framework realizing the pinned subspace-incidence 
system {H,X) is a triple {H,X,p), where p is a realization of i?, such that for 
all pins Xk € X, Xk is contained in {p(efe)), the subspace spanned by the set of 
points realizing the vertices of the hyperedge corresponding to Xk- 

We may write m{xk) or simply mk in substitute of m{ek), where Xk is the 
pin associated with the hyperedge Cfc. 

Since we only care about incidence relations, we projectivize the Euclidean 
space to treat the pinned subspace-incidence system in the real projective 
space P‘^“^(]R), and use the same notation for the pins and hypergraph realiza¬ 
tion when the meaning is clear from the context. 

Figure [lb] gives an example of a pinned subspace-incidence framework in the 
projective space with d = 4, where each pin is the subspace spanned 

by the set of cross-denoted points on the corresponding hyperedge. 

Definition 2. A pinned subspace-incidence system {H,X) is independent if 
none of the algebraic constraints is in the ideal generated by others, which 
generically implies the existence of a realization. It is rigid if there exist at 
most finitely many realizations. It is minimally rigid if it is both rigid and 
independent. It is globally rigid if there exits at most one realization. 


4 Algebraic Representation and Linearization 

In the following, we use A[R, C] to denote a submatrix of a matrix A, where 
R and C are respectively index sets of the rows and columns contained in the 
submatrix. In addition, A[R, • ] represents the submatrix containing row set R 
and all columns, and A[- ,C\ represents the submatrix containing column set 
C and all rows. 

4.1 Algebraic Representation 

A pin Xk associated with a hyperedge Ck = {nf,f 2 , •.. of cardinality Sk 
is constrained to be contained in the subspace {p{ek)) spanned by the point 
set {pi,P 2 , ■ ■ ■ iPsj,}. As Xk is a subspace of dimension mk — 1 in P^“^(]R), we 
can pick a set of points {x\,x\^... , 2 ;^^.} spanning Xk from Xk- Now the 
constraint is equivalent to requiring each such point xf to lie on {p(e)), for 
1 < ^ < TOfc. We call each such point xf a multipin as it acts like a pin with 

m{xi) = 1 . 

Using homogeneous coordinates pf = [ pf^i ■ ■ ■ Pi,d-i ] xf = 

[ x^ x ^2 ■ ■ ■ ^fd -1 ]j write this incidence constraint for each point xi 
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by letting all the Sk x Sk minors of the Sfc x (d — 1) matrix 

k 

Pi ~ 
p^-xf 

pI - . 

be zero. There are minors, giving equations. Note that any d — Sk 
of these equations are independent and span the rest. So we can write 

the incidence constraint as {d — Sk) independent equations: 

det(i;f[-,C(t)]) =0, l<t<d-Sfe (1) 

where C{t) denote the following index sets of columns in E: 

C{t) = {1, 2,... Sfe - 1} U {sfe - 1 + t} 1 < t < d - s 

In other words, C{t) contains the first Sfe — 1 columns together with Column 
Sk 1 “t“ t. 

Now the incidence constraint for the pin Xk is represented as mk{d — Sk) 
equations for all the mk multipins {x\^X 2 , ■ ■ Consequently, the pinned 

subspace-incidence problem reduces to solving a system of ~ ^k) 

equations, each of form ©■ We denote this algebraic system by {H,X){p) = 0. 

4.2 Genericity 

We are interested in characterizing minimal rigidity of pinned subspace-incidence 
systems. However, checking independence relative to the ideal generated by the 
variety is computationally hard and best known algorithms, such as computing 
Grobner basis, are exponential in time and space m- However, the algebraic 
system can be linearized at generic or regular (non-singular) points, whereby 
the independence and rigidity of the algebraic system {H,X){p) = 0 reduces to 
linear independence and maximal rank at generic frameworks. 

In algebraic geometry, a property being generic intuitively means that the 
property holds on the open dense complement of an (real) algebraic variety. 
Formally, 

Definition 3. A pinned subspace-incidence system {H,X) is generic w.r.t. a 
property Q if and only if there exists a neighborhood N{X) such that for all 
systems {H,X') with X' G N(X), (H,X') satisfies Q if and only if {H,X) 
satisfies Q. 

Similarly, a framework {H,X,p) is generic w.r.t. a property Q if and only 
if there exists a neighborhood N(p) such that for all frameworks (7J, A, q) with 
q G N{p), {H,X,q) satisfies Q if and only if {H,X,p) satisfies Q. 

Furthermore we can define generic properties in terms of the underlying 
weighted hypergraph. 
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Definition 4. A property Q of pinned subspace-incidence systems is generic 
(i.e, becomes a property of the underlying weighted hypergraph alone) if for 
any weighted hypergraph H = (V,E,m), either all generic (w.r.t. Q) systems 
{H,X) satisfies Q, or all generic (w.r.t. Q) systems {H,X) do not satisfy Q. 

Once an appropriate notion of genericity is defined, we can treat Q as a prop¬ 
erty of a hypergraph. The primary activity of the area of combinatorial rigidity 
is to give purely combinatorial characterizations of such generic properties Q. 
In the process of drawing such combinatorial characterizations, the notion of 
genericity may have to be further restricted by so-called pure conditions that 
are necessary for the combinatorial characterization to go through (we will see 
this in the proof of Theorem fTTl) . 

4.3 Linearization as Rigidity Matrix 

Next we follow the approach taken by traditional combinatorial rigidity theory 
mm to show that rigidity and independence (based on nonlinear polynomials) of 
pinned subspace-incidence systems are generically properties of the underlying 
weighted hypergraph H, and can furthermore be captured by linear conditions 
in an infinitesimal setting. Specifically, LemmalHlshows that rigidity of a pinned 
subspace-incidence system is equivalent to existence of a full rank rigidity matrix^ 
obtained by taking the Jacobian of the algebraic system {H,X){p) at a regular 
point. 

A rigidity matrix of a framework (i/. A, p) is a matrix whose kernel is the 
infinitesimal motions (flexes) of {H,X,p). A framework is infinitesimally rigid 
if the space of infinitesimal motion is trivial, i.e. the rigidity matrix has full 
rank. To define a rigidity matrix for a pinned subspace-incidence framework 
{H, X,p), we take the Jacobian of the algebraic system (iJ, X){p) = 0 by taking 
partial derivatives with respect to the coordinates of pi’s. In the Jacobian, each 
vertex Vi has d — 1 corresponding columns, and each pin Xk associated with the 
hyperedge Cfc = {nf, V 2 , ■ ■ ■, has mk{d—Sk) corresponding rows, where each 
equation det (Af [ • , C{t)]) = 0 ([T]), i.e. Equation t of the multipin xf, gives the 
following row (the columns corresponding to vertices not in Ck are all zero): 


9 det {Ef[ • 

,C{t)]) 

a det {E^ [ • 

,c{t)]) 

a det (Af)- 

,c{t)]) 

n 0 

dp\^ 

1 

dp\ 

1 • ■ 

2 

dpU- 

-1 

, u, u 

a det (Af[- 


a det [Ef [ • 

,Cit)]) 

a det (Af [ • 


n 0 

dp2, 

1 

dpi 

5 • ■ 

2 

dpid- 

-1 

, u, u 

a det (Af[- 


a det [Ef [ • 

,cm 

a det (Af)- 


0 

dPsk. 

,1 

dpi. 

•> • ■ 

,2 

9pL<i 

-1 

, u,. 


(2) 
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Let be the matrix whose rows are coordinates oi Pi,P 2 , ■ ■ ■ ,P^^' 


Pl.l 

Pi,2 ■ 

• Pl,d-1 

P2.1 

P2,2 

■ P2,d-1 

. Psk,l 


Psh,d-1 


Let V^^ be the V^[- ,C'(t)], i.e. the Sk x Sk submatrix of containing only 
columns in C{t). Let be the matrix obtained from by replacing the 
column corresponding to Coordinate j with the all-ones vector (1,1,...,!) for 
j £ C{t), and the zero matrix for j ^ C{t). Let be the determinant of 
Let xi = (note that = !)• © can be rewritten in 

the following simplified form: 


0 , . 

.., 0 , 0 , Dl.bWDl^hW ■ 

Dk yk,l 

0 , 




.., 0 , 0 , Dl^bl'\ ■ 


0 , . 




.., 0 , 0 , . 

■ 5 ^t,d-l^Sk ’ 

0 , . 

.,0 

(3) 


Each vertex has the entries jb^’\ 1 < j < d — 1 in its d — 1 columns, 
among which exactly Sk entries with j £ C(t), i.e. the first S/t — 1 columns 
together with Column — 1 +1, are generically non-zero. Note that the terms 
are equal for all t, so we may just use to denote it. 

For each 1 < t < d — Sk, there are ruk rows as where each multipin xf 
corresponds to the row r*; for 1 < I < ruk- These ruk rows have exactly the 
same row pattern except for different 6^’*’s: 


Vl.l 

Vl.2 

VlA-l 



■ * ,d—l 

DlA'" 

D^bl'" ■ 

jjk Lfe.l 


jjk Lfc.l 
^t,2^Sk 

jjk uk,l 

^t,d-l^sk 



T^k ^k,2 


-^t,2^Sk 

jjk uk,2 

^t,d-l^Sk 

nk lk,mk 

nk lk,mk 

k-'t,2^1 

nk ik,'mk 

■ ■ • 

Dk uk,mk 
^t,2^Sk 

j~)k ]-,k,mk 

■ ■ ^t,d-l°Sk 


Example 1. For d = 4, consider a pin x with m{x) = 2 associated with the 
hyperedge e = {vi,V 2 \- The pin has the following m(x) ■ (d — s(e)) = 4 rows in 
the simplified Jacobian (the index k is omitted): 




Vip 

Vl,2 

fl,3 

V2,l 

V2,2 

V2,3 

t = 

1,1=1 

Diib\ 

Dbl 

0 

Di^ib\ 

Dbl 

0 

t = 

1,1 = 2 

... 

Dbl 

0 


Dbl 

0 

t = 

2,1 = 1 

7^2,1^1 

0 

Dbl . 

D2,ib\ 

0 

Dbl 

t = 

2,1 = 2 

- ... D2,lbl 

0 

Dbl ■ 

■ D2,lbl 

0 

Dbl 
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We define the rigidity matrix M {H, X, p) or simply M (p) for a pinned 
subspace-incidence framework {H, X,p) to be the simplified Jacobian matrix ob¬ 
tained above, where each row has form (|31). It is a matrix of size mk{d — Sk) 
by n{d- 1 ). 

Definition 5. A pinned subspace-incidence framework {H,X,p) and the cor¬ 
responding system {H,X) is generic Up and X are regular/non-singular points 
with respect to the algebraic system {H,X){p) — 0. 

We use M{H) of simply M to denote the generic rigidity matrix for a 
weighted hypergraph H. Note that the rank of M cannot be less than the 
rank of M (p) for any specific realization p. 

Lemma 6. Generic infinitesimal rigidity of a pinned subspace-incidence frame¬ 
work {H,X,p) is equivalent to generic rigidity of the system {H,X). 

The proof of Lemma [5] follows the traditional combinatorial rigidity ap¬ 
proach [T] and is given in the Appendix. 


5 Combinatorial rigidity characterization 

5.1 Required hypergraph properties 

This section introduces pure hypergraph properties and definitions that will be 
used in stating and proving our main theorem. 

Definition 7. A hypergraph H = {V,E) is (/c, 0)-sparse if for any V' C V, 
the induced subgraph H' = {V',E') satisfies \E'\ < k\V'\. A hypergraph El is 
(fc,0)-tight if H is (fc, 0)-sparse and \E\ = k\V\. 

This is a special case of the (fc, Z)-sparsity condition that was formally studied 
widely in the geometric constraint solving and combinatorial rigidity literature 
before it was given a name in m- A relevant concept from graph matroids is 
map-graph, defined as follows. 

Definition 8. An orientation of a hypergraph is given by identifying as the 
tail of each edge one of its endpoints. The out-degree of a vertex is the number 
of edges which identify it as the tail and connect v to V — v. A map-graph is a 
hypergraph that admits an orientation such that the out degree of every vertex 
is exactly one. 

The following lemma from [15] follows Tutte-Nash Williams iniiii] to give 
a useful characterization of (fc, 0 )-tight graphs in terms of maps. 

Lemma 9. A hypergraph El is composed of k edge-disjoint map-graphs if and 
only if H is {k,0)-tight. 

Our characterization of rigidity of a weighted hypergraph H is based on 
map-decomposition of a multi-hypergraph H obtained from H. 
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Definition 10. Given a weighted hypergraph H = {V,E,m), the associated 
multi-hypergraph H = {V, E) is obtained by replacing each hyperedge e/c in E 
with a set E’^ of mk{d — Sk) copies of multi-hyperedges. 

A labeling of a multi-hypergraph H gives a one-to-one correspondence be¬ 
tween E^ and the set of mk{d — Sk) rows for the hyperedge in the rigidity 
matrix M, where the multi-hyperedge corresponding to the row r^ ^ is labeled 



5.2 Characterizing rigidity 

In this section, we apply m to give combinatorial characterization for minimal 
rigidity of pinned subspace-incidence systems. 

Theorem 11 (Main Theorem). A pinned subspace-incidence system is generi- 
cally minimally rigid if and only if: 

I C'l 

1. The underlying weighted hypergraph H = {V,E,m) satisfies rnk{d — 

Sk) = (d- 1)|F|, and Y)ek&E' '^k{d - Sk) < {d - for every vertex 

induced subgraph H' = {V',E'). In other words, the associated multi- 
hypergraph H = (V, E) has a decomposition into {d — 1) maps. 

2. There exists a labeling of H compatible with the map-decomposition (de¬ 
fined later) such that in each set E^ of multi-hyperedges, 

(a) two multi-hyperedges e^^ and e^^ with li — I 2 are not contained in 
the same map in the map-decomposition, 

(b) two multi-hyperedges e^^ and e^^ with ti = t 2 do not have the same 
vertex as tail in the map-decomposition. 

To prove Theorem 1111 we apply Laplace expansion to the determinant of 
the rigidity matrix M, which corresponds to decomposing the (d — l,0)-tight 
multi-hypergraph id as a union of d — 1 maps. We then prove det(M) is not 
identically zero by showing that the minors corresponding to each map are not 
identically zero, as long as a certain polynomial called pure condition is avoided 
by the framework. The pure condition characterizes the non-genericity that the 
framework has to avoid in order for the combinatorial characterization to go 
through: see Example ID 

A Laplace expansion rewrites the determinant of the rigidity matrix M as a 
sum of products of determinants (brackets) representing each of the coordinates 
taken separately. In order to see the relationship between the Laplace expansion 
and the map-decomposition, we first group the columns of M into d — 1 column 
groups Cj according to the coordinates, where columns for the first coordinate 
of each vertex belong to Ci, columns for the second coordinate of each vertex 
belong to C 2 , etc. 

Example 2. For d = 4, consider a pin x with m{x) = 2 associated with the 
hyperedge e = vi,V 2 . The regrouped rigidity matrix has d — 1 = 3 column 
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groups, where the pin x has the following 4 rows (the index k is omitted): 





Vip 

V 2.1 

■Cl ,2 

V2,2 

Vl,3 

V2,3 


t = 

1 

1 = 1 

■ ... Di^ib] . 

■ Di^ib\ ... 

... Dbl . 

. Dbl 



- 

t = 

1 

/ = 2 

■ ■ ■ Di,ibi ■ 

■ Di ib2 ... 

... Dbl ■ 

. Dbl ■■■ 




t = 

2 

1 = 1 

■ ■ ■ D2pb\ . 

■ D2pb2 ... 



... Db\ 

. Dbl 


t = 

2 

1 = 2 

. ... D2,lbl . 

■ L*2,1^2 ■ ■ • 



... Dbl ■ 

. Dbl ■ 



We have the following observation on the pattern of the regrouped rigidity 
matrix. 

Observation 1. In the rigidity matrix M with columns grouped into column 
groups, a hyperedge Cfc has mk{d — Sk) rows, each associated with a multi¬ 
hyperedge of e^ in H. In a column group j where j < Sfe — 1, each row associated 
with Ck contains Sk nonzero entries at the columns corresponding to vertices of 
Cfe. In a Column group j where j = Sfc — 1 + t > Sfc, a row assocated with a 
multi-hyperedge ; of e^ is all zero if r ^ t; the remaining mk rows contains 
Sk nonzero entries at the columns corresponding to vertices of Ck- 

A labeling of H compatible with a given map-decomposition can be obtained 
as following. We start from the last column group of M and associate each col¬ 
umn group j with a map in the map-decomposition. For each multi-hyperedge 
of the map that is a copy of the hyperedge e^, we pick a row r^j that is not 
all zero in Column groups j and label the multi-hyperedge as e^ j. By Observa¬ 
tion [Jl this is always possible if each map contains at most mk multi-hyperedges 
of the same hyperedge Ck, which must be true if there exists any labeling of H 
satisfying Condition 2(a) of Theorem fTTl 

In the Laplace expansion 


det(M) = 

(7 


±]^detM[i?' 



(4) 


the sum is taken over all partitions a of the rows into d — \ subsets Rf,R 2 , 
..., i?J, ..., Rfi_i, each of size \V\. In other words, each summation term 
of ([3]) contains |I4| rows Rf from each column group Cj. Observe that for any 
submatrix M[Rj , Cj], each row has a common coefficient Dtj, so 


det{M[R'^,Cj]) 


[ n DlAdet{M'[R^,C,]) 


where each row of M'[Rj ,Cj] is either all zero, or of the pattern 


[0,...,0,b1-\b^2\0> 


0 01 


with non-zero entries only at the Sk indices corresponding to vf € Cfc. 
For a fixed a, we refer to a submatrix M[RJ,Cj] simply as Mj. 


(5) 
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Figure 2: (a) A minimally rigid pinned subspace-incidence system in d = 3. 
(b) A map-decomposition of the multi-hypergraph of the system in (a), where 
multi-hyperedges with different colors are in different maps, and the tail vertex 
of each multi-hyperedge is pointed to by an arrow. 


Example 3. Figure [2^ shows a pinned subspace-incidence system in d = 3 with 
4 vertices and 5 hyperedges, where ei = {vi},e 2 = {r' 2 },e 3 = {^ 1 ,^ 3},64 = 
{v2, ^^4}, 65 = {^3, U4}, and irik = 1 for all 1 < fc < 5 except that = 2. Figure 
I2bl gives a map-decomposition of the multi-hypergraph H of (a). The labeling of 
multi-hyperedges is given in the regrouped rigidity matrix where the colored 
rows inside the column groups constitute the submatrices MJ in the summation 
term of the Laplace decomposition corresponding to the map-decomposition. 
The system is generically minimally rigid, as the map-decomposition and label¬ 
ing of H satisfies the conditions of Theorem [TT] 


e 

e 

e 

e 

e 

e 

e 

e 


1 

1.1 

1 

2.1 

2 

1.1 

2 

2.1 

3 

1.1 

4 

1,1 

5 

1,1 

5 

1.2 



( 6 ) 

Main Theorem. First we show the only if direction. For a generically minimally 
rigid pinned subspace-incidence framework, the rigidity matrix M is generically 
full rank, so there exists at least one summation term a in where each 
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submatrix Mj is generically full rank. As the submatrices don’t have any over¬ 
lapping rows with each other, we can perform row elimination on M to obtain 
a matrix N with the same rank, where all submatrices Mj are simultaneously 
converted to a permuted reduced row echelon form Nj , where each row in Nj has 
exactly one non-zero entry at a unique Column i, In other words, all iVj’s 
can be converted simultaneously to reduced row echelon form by multiplying a 
permutation matrix on the left of N. Now we can obtain a map-decomposition 
of H by letting each map j contain multi-hyperedges corresponding to rows of 
the submatrix Nj, and assigning each multi-hyperedge in map j the vertex i 
corresponding to the non-zero entry jdf in the associated row in Nj as tail. In 
addition, such a map-decomposition must satisfy Condition 2 in Theorem 1111 

Condition 2(a): assume two multi-hyperedges ; and ; are in the same 
map j, i.e. the rows corresponding to these two edges are included in the same 
submatrix Mj. If j > s — 1, one of these rows must be all-zero in Mj by 
Observation[T] contradicting the condition that Mj is full rank. If j < s — 1, both 
of these rows in Mj will be a multiple of the same row vector ([5]), contradicting 
the condition that Mj is full rank. 

Condition 2(b): note that the rows in M corresponding to multi-hyperedges 
and have the exactly the same pattern except for different values of 
6’s. If e^i^ has vertex i as tail, after the row elimination, the column containing 
will become the only non-zero entry of column i for iVjj, while the column 
containing will become zero in all column groups, thus i cannot be assigned 
as tail for 

Next we show the if direction, that the conditions of Theorem HD imply 
infinitesimal rigidity. 

Given labeled multi-hypergraph H with a map-decomposition satisfying the 
conditions of Theorem m we can obtain summation term a in the Laplace 
decomposition (|4]) according to the labeling of H, where each submatrix Mj 
contain all rows corresponding to the map associated with Column group j. 

We first show that each submatrix Mj is generically full rank. According 
to the definition of a map-graph, the function t : E ^ V assigning a tail 
vertex to each multi-hyperedge is a one-to-one correspondence. We perform 
symbolic row elimination of the matrix M to simultaneously convert each Mj 
to its permuted reduced row echelon form Nj, where for each row of Nj, all 
entries are zero except for the entry corresponding to the vertex T{e^i), 
which is a polynomial in b^'^'s in the submatrix Mj. Since Mj cannot contain 
two rows with the same k and I by Condition 2(a), the in different rows of a 
same map are independent of each other, 0 under a generic specialization 

of b^'K Since each row of Nj has exactly one nonzero entry and the nonzero 
entries from different rows are on different columns, the \V\ x \V\ matrix Nj is 
clearly full rank. Thus Mj must also be generically full rank. 
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(S') (b) 

Figure 3: (a) A pinned subspace-incidence system in d = 4. (b) A map- 

decomposition of the multi-hypergraph of the system in (a), where multi¬ 
hyperedges with different colors are in different maps, and the tail vertex of 
each multi-hyperedge is pointed to by an arrow. 


We conclude that 

det(M)=^ [±n(( n (7) 

where the sum is taken over all a corresponding to a map-decomposition of 
H. Generically, the summation terms of the sum © do not cancel with each 
other, since det(M'[i?J, C^]) are independent of the multi-linear coefficients 
Hr'' PR'' rows of M are independent by Condition 2(b). This 

implies that M is generically full rank. 

The polynomial ([7]) gives the pure condition for genericity. In particular, 
when there is a subgraph {V, E') with \V'\ < d and ^ 

condition vanishes and the system won’t be minimally rigid: see Example 0] 

□ 


5.2.1 Pure condition 

The pure condition ([7]) obtained in the proof of Theorem [TT] characterizes the 
badly behaved cases that break the combinatorial characterization of infinites¬ 
imal rigidity. However, the geometric meaning of the pure condition is not 
completely clear. One particular condition not captured by Theorem [TT] but 
enforced by the pure condition is that there cannot exists a subgraph {V',E') 
of H with \V'\ < d such that Xe^gS'^^ ^ otherwise simple counterex¬ 
amples can be constructed to the characterization of the main theorem. An 
immediately consequence is that for any hyperedge e^, the dimension ruk of its 
associated pin must be less than or equal to its cardinality Sk- 

Example 4. Figure [3^ shows a pinned subspace-incidence system in d = 4 with 
4 vertices and 5 hyperedges, where rrik is 2 for fc = 5 and is 1 otherwise. A map- 


13 








decomposition of the multi-hypergraph H of the system is given in Figure I3b[ 
and we can easily find a labeling of H satisfying conditions in Theorem 1111 
However, the system is not minimally rigid, as generically the pin Xi will not fall 
on the plane spanned by pins X 4 and X 5 . Note that the sub-hypergraph (H', E') 
spanned by vertices vi,V 2 ,V 4 violates the pure condition as X^e^G-E' = itii + 
m 4 -I- TO 5 = 4 > \V'\ = 3. 

6 Conclusion 

In this paper, we studied the pinned subspace-incidence, a class of incidence 
geometric constraint system with applications in dictionary learning and bio¬ 
material modeling. We extend our results in m and obtain a combinatorial 
characterization of minimal rigidity for general pinned subspace-incidence sys¬ 
tems with non-uniform underlying hypergraphs and pins being subspaces with 
arbitrary dimensions. 

As future work, we plan to extend the underlying group of the pinned 
subspace-constraint system, i.e. consider two frameworks to be congruent if the 
point realization and pin set of one can be obtained from the other under the 
action of a certain group, for example the projective group. Another possible 
direction is to apply Cayley factorization to find geometric interpretations 
of the pure conditions. 
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A Proof of Lemma [6] 


Proof. First we show that if a framework {H, X,p) is regular, infinitesimal rigid¬ 
ity implies rigidity. Consider the polynomial system (H,X){p) of equations. 
The Implicit Function Theorem states that there exists a function such that 
p = g{X) on some open interval, if and only if the rigidity matrix M has full 
rank. Therefore, if the framework is infinitesimally rigid, the solutions to the 
algebraic system are isolated points (otherwise g could not be explicit). Since 
the algebraic system contains finitely many components, there are only finitely 
many such solution and each solution is a 0 dimensional point. This implies 
that the total number of solutions is finite, which is the definition of rigidity. 

To show that generic rigidity implies generic infinitesimal rigidity, we take 
the contrapositive: if a generic framework is not infinitesimally rigid, we show 
that there is a finite flex. If {H,X,p) is not infinitesimally rigid, then the rank 
r of the rigidity matrix M is less than {d— 1)|I^|. Let E* be a set of edges in 
H such that |i?*| = r and the corresponding rows in M are all independent. 
In M[E*, •], we can find r independent columns. Let p* be the components 
of p corresponding to those r independent columns and p*'^ be the remaining 
components. The r-by-r submatrix M[E*,p*], made up of the corresponding 
independent rows and columns, is invertible. Then, by the Implicit Function 
Theorem, in a neighborhood of p there exists a continuous and differentiable 
function g such that p* = g{p*'^). This identifies p', whose components are 
p* and the level set of g corresponding to p*, such that {H,X){p') = 0. The 
level set defines the finite flexing of the framework. Therefore the system is not 
rigid. □ 
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